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Abstract 

<N 

In this paper, we consider zero range process with an initial condition which is equiv- 
{N| alent to step initial condition in total asymmetric simple exclusion process (TASEP) 

as described in a paper by Rakos, A. and Schiitz p2]by using techniques developed by 
Borodin, Ferrari, and Sasamoto [2]. The solution for the transition probability of total 
asymmetric simple exclusion process for particles with different hopping rates was first 
worked out by Schiitz and Rakos [TT] for the case when p = 1 or q = 1. The formula 
was later applied to analyze two speed TASEP [2] with alternating initial condition. 
Here we will investigate the two speed TASEP case with step initial condition. 

J: 

1 Introduction 
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The asymmetric simple exclusion process (ASEP) was first introduced in 1970 by Frank 
Spitzer [TJ] . It has a wide range of applications [9 J [14] and it is now known to belong to the 
tS- Kardar-Parisi-Zhang (KPZ) universality class [7]. 

ASEP is a model for particles interacting on a lattice. Each particle is equipped with an alarm 
clock with exponential distribution of hopping rate vf The particle jumps when the clock 
rings. It has a probability p jumps to the right and q = 1 — p jumps to the left provided the 
neighboring site is empty. When p — 1, we call it total asymmetric simple exclusion process 
(TASEP). Kurt Johansson[5] was the first to show that when N particles with equal hopping 
rates are initially assigned in the negative region called step initial condition, meaning the 
^ position of the ith particle Xi(t) = — i + 1 when t — 0, the following formula connecting 

TASEP with random matrix theory holds: 



P(x N (t)>t)= / ... / l\t^ N - 1 e-^T[(t l -t 3 ) 2 dt 1 ---dt N . 

J ° J ° i=l i<j 



Here, the right hand side of (1.1.1) with the Laguerre weight function (w(U) = t^ +N ~ 1 e^ ti ), 
is equal to the probability that the largest eigenvalue of a random matrix AA* is < t in 
Laguerre ensemble where A is a N x (x — 1 + 2iV) matrix of complex Gaussian random 
variables with mean zero and variance 1/2. 

Later on Gunter M. Schiitz and A. Rakos [TT] and Tomohiro Sasamoto and Taro Nago [13] 
also derived the same formula above using Bethe Ansatz technique. Recently Craig Tracy 
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and Harold Widom [16] have also shown that Bethe Ansatz technique can be applied to 
ASEP process. They also extended Johansson result to ASEP process. 

In this paper, we will use the Fredholm determent formula for P(xi(t), ...,Xx(t);yi, ....,yjy) 
from a paper by Borodin and Ferrari [TJ. Then we will work out the bi-orthornormal poly- 
nomials for the correlation kernel of P(xi(t), Xjv(t); yi, yjsi)- We will see that when 
the particles starts with step initial condition while the first M leading particles moving 
with rate a G (0, 1), slower than the rest of the particles which have unit speed, our process 
diagram, as shown in section 2, has a frozen regiorj^] and it is different from the case with 
periodic initial condition [2J. Also with step initial condition, we won't have a so called shock 
line as the case in the paper by Borodin, Ferrari, and Sasamoto [2]. We will also discuss how 
we can apply the formula in section 3 to study zero range process as certain parameters goes 
to infinite as what Rakos, A. and Schiitz discussed in their paper 



2 Main Results 

We will consider TASEP With step initial condition Xj(0) = y^ = M — i and letting the first 
M particles have rate a and the rest of particles with unit speed(see figure below for an 
example with M = 3). 
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In this section, we will study the asymptotics of P(x n (t) > x). The main results to be proved 
in this section is the following process diagram which gives the relationship between n/t and 
a. Here n and t are both parameters in x n (t) and a is the speed of the right most particle. 
In the diagram is the Airy2 process, Adbm is the Dyson's Brownian motion process and 
the curved line (1 — a) 2 in the diagram represents the Adbm-^2 process. The definition of 
these process will be given in appendix. It should be noted here our diagram has a frozen 
region and it is different from the case with periodic initial condition[2]. 



1 Frozen region means particles sitting in this region are standing still since their neighboring sites are all 
occupied. 
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n/t 




Cruve ine is the A_DBM to 2 process 
Horizontal line is the Frozen region 
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A by product of our result is to use our results to study zero range process which will be 
explained in the following paragraphs. 

Zero range process is a process similar to TASEP but here our particles do not exclude 
each other and all particles' hopping rates depend on their which site they are sitting on. 
In our case we will assume every particle has rate 1 hopping to the left from site 1 to site 
L — 1 and there is an injection of particles to site L — 1 from the right with rate a. We will 
call this kind of zero range process totally asymmetric zero range process. It is easy to see 
the current distribution across the (L - l) m bond of our total zero range process (TAZRP), 
PL_i(x,t), is equal to the (L)-particle TASEP case F(xi,+i(t) > x). When we map the L 
sites zero range process to (L + 1) particles TASEP, we have Vi — 1 for all i but v l — ex. 
This approach of using Bethe ansatz solution to study zero range process was also used by 
Rakos and Schutz[T2]. In this paper, we take the limit L = vt then we have 




(2) 



The picture of mapping zero range process to TASEP is given below. 
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Now we will discuss our model and our main results. If particles starts from Z_ then the 
particles density P(there is a particle at x at time t):= u(x,t) is given as below [8] 

1, x < —t, 

u(x,t) = { l/2-x/t/2, x G [— t, (2a — l)t], (3) 
1 — a, x > (2a — l)t. 

while the initial condition is given as the following 

u(x,0):=\)> ( 4 ) 
1 — a, x > 0. 
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u(x,0) 



Ll(X.t) 
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(1-a) 



The particle at position at is the right most particle with rate a. 

By Burke's Theorem jl], when there are M particles with rate a placed on N initially, our 
initial condition is identical to step initial on Z_ and Bernoulli with density (1 — a) on Z+. 
Mathematically, the continuous version of u(x,t), let's say p(£, r) with p(x/t,l) = u(x,t), 
satisfies the Burger's equation 

d T p + dt(p(l-p)) = 0. (5) 
Since u(x,t) = 1 — a when x > (2a — l)t is a constant and there are (1 — a) 2 t particles 



moving around speed a , we know if n E (M. 
to Xuif) — oct which is n — M = f* M u(x, t)dx 



a) i t), n = area under u(x,t) from x n (t) 



n — M = (at — x n )(l — a) 

^x n = at- (n- M)/(l- a). (6) 

If we apply the above steps to the other regions and let n/t = v, then we have 

E(rt t f \\ [ oc-u/(l-a), v G (0, (1 — a) 2 ), 
ton HX[ " m := 1-2^, ue((l-ani), (7) 



t— >OD 



—v, v G (1, oo). 



Please be noted that we will always assume < a < 1. 

Since the particles get slowed down when v e [0, (1 — a) 2 ], we will look at fluctuations on a 
\Jt scale in this region . For v G ((1 — a) 2 , (1)), we will look at fluctuations on a i 1 / 3 scale, 
for v > 1, there is no fluctuation. 



2 Here we use ¥,(x n (t)) = x n (t) as the notation for the average position of the particle x n (t) at time t. 
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We now have the rescaled process given as 

f Xfwr{<rf :^ )/(1 " a)} , - e (0, (1 - «) 2 ), a» = a(l - u/(l - af) 
X t (u):=l ^zftzgVM) , z,G((l-«) 2 ),l), ( 8 ) 

[ o, i/ e (l, oo). 

We now let 7r(0) be a function from R — > R such that | vr' j < 1. Define £j = (tt(^) + 9j)T , 
rii = {n{0i) — 6i)T + M , a = n — 9, and u = it + 9 then 



W) (^M^-^) ; ^ G (Q) (1 _ q)2)j ^ = ^ _ a/(1 _ a)2 

o, i/ e (l, oo). 

(9) 

Below is the process diagram which we will prove in this chapter . 
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Process Diagram: 




v 



In section 4, we will sketch the proof of convergence of the rescaled process X? in each region 
in the sense of finite dimensional distributions. 



3 Kernels 

3.1 Overview 

In this section, we briefly discuss about the correlation kernel for the joint probability dis- 
tribution of the position for arbitrary number of particles. 

Here we will write the joint probability distribution of positions Xi(t) as a Fredholm de- 
terminant expression. Before doing that, we first consider a set of non-decreasing positive 
numbers {v\, ...,v n } such that Vi < Vj if i < j. Let {u\ < «2 < ••• < u u } be their different 
values with a k being the multiplicity of u k . Also we define a set called space like set 

S = {(n k , t k ),k = 1, m\n k > n k+1 , t k < t k+1 }. 

Define 

V n = {x l u%\l < k < z/,0 < I < a k - 1} (10) 
as a space of functions in x. 



Theorem 3.1. 1[T]: 

Consider particles such that Xj(0) = Take a sequence of particles and times in the space 
like set S. The joint probability distribution of x nk (t k ) is given by the following Fredholm 



7 



determinant expression 

P{r^ =x {x nk {t k ) > a k }) = det(l - Xa-K'Xo){j3({(ni,t 1 ) I ...,(n m ,t m )}xZ)} (11) 

where Xa(( n fc, £&), x) = h x <a k }- Here if is a kernel on l 2 (S x Z) and it is defined as the 
following 



K((n u h), Xl ; (n 2) t 2 ),x 2 ) = -<P^^\x 1} x 2 ) + £ Ki*M*%-kM ( 12 ) 

fe=i 

where 

1 / dw e tw Y\k=i( w ~ v k) 

J '"' ! 27Ti 



and 

.((n^.^W X2 ) = ±l ^ e(tl " t2)W I* (14) 

2m J T w w x+ni ' X2 - U2 (w - l) n2 - ni 

where I5 is the identity function on the space like space S , v — {v ni+ x, v n2 ) and T ^ 
denotes any positive oriented simple closed curve that includes and the points in the set v. 

The functions {®n-j}i<j<n are described by the following conditions: 

<^t j ,K'- k >=S j , k , l<j,k,<n, (15) 
and span{$™^(x), 1 < j < n} = V n . 

The orthornormal polynomials for step initial condition i/j = M — j are given at next section. 



3 Another way to see Fredholm determinant is 

det(l - Xa-^Xo)/ 2 (({(n 1) ii),...,(n m ,i TO )})xZ) 

= Y~^~ Y Y "' Y detK({n ik ,t ik ,x k ); (« 4j ,ii J ,a;j))|i<i,fc<„ 

n>0 i\ ,. . . 7 in — 1 xi >ai x n ~>a n 

= Y ~ Y / <%•••/ dy n detii'((ni Ji ,ti fc! [yfc]);(n^,t i;!! [y i ]))|i<j !fe <„ 

n>0 i 1 ,-,in=l ,/{l ' 1 - 0l l J{y n >a n } 

if [j/j] = Lowest integer greater or equal to yi. 



3.2 Bi-orthornormal polynomials 



In this section we verify the orthogonal condition for our kernel with step initial condition 
(see picture below). 
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Orange particles have hopping rate a and black particles have hopping rate 1. 

The main difficulty to apply Theorem 3.1.3 is to find the corresponding bi-orthogonal poly- 
nomials for its kernel. It was pointed in [2] that bi-orthogonally polynomials in [2] for the 
case of alternating initial condition can be obtained by Gram-Schmidt orthogonalization 
procedure. However, we do not do it here. We obtained our bi-orthogonal polynomials by 
modifying the bi-orthogonal polynomials for the case of alternating initial condition. They 
are determined once the bi-orthogonal condition is meet. 

We will show that the &™-j( x ) functions which we will define later satisfies the following 
relationship 

< K'-p K'-k >■= E = *J* 1 < 3, *, < n. (16) 

and _ k are called bi-orthogonal if they all satisfy the above relationship. The above 
<, > map can be viewed as an inner product function from V n x V n — > Z. Using the Cauchy 
integral formula, one can see the $^ -(ar) functions span V n . Thus, together with the following 
theorem, wm satisfy the condition for Theorem 3.1.1. 

From now on, every contour integral's [^j differential dz is identified as 

In previous section, we know our kernel is expressed as a pair of bi-orthogonal polynomials 
$ and our polynomials is constructed based on the polynomials in [2]. We will follow 
their notations as much as possible. With step initial condition yi = M — i and letting the 
first M particles have rate a and the rest of particles with unit speed(see figure above for an 
example with M = 3), the orthonormal polynomials for the correlation kernel in previous 
section are given as the following which is slight different from the case of alternating initial 
condition [2]. 

Theorem 3.2.1 Bi-orthonormal polynomials for the Kernel K: 

From now on we use the notation r pi) ... )Pj to represent a positive oriented contour integral 
enclosed the points: p%, ...,p n . 



instead of writing § we write § dz. 
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(a) For n< M, 



1_ f dw (w - a) n - J 
" " ' 2:;i % w w x+n ~ M 



$ n '* .( x ) = — <t dv {i+v y +n M e -t(v+i) 

2tt* /r^ - (« - 1))^ +1 



(b) For n > M and j > M + 1, 



n JK ' 2ni J Toi w w x+n ~ M 

n jK ' 2-Kl J To V n-J+1 



(c) For n > M and j < M, 



n ~ 3y ' 2-Kl J Foa W w x+n-M 

nt i r i r i (i + z ) x+n - M 

V ' llJ ~ (2niy % a _ x dV (v -(a- l))M-J + i % Qv dz z _ „ z n-M e t iz+ i) 6 
Here, 

( ) 2z + A ( Z + A ) M ~ 3 
9[z ' v) - Z + V + A (v + A) M ~^ 

and A » \z\,z G To,*,. 
Proof: 

First we recall by definition 

x& 

and since ■(#) = ii x < M — n — lwe have 

x£Z x>M-N 

The proof is by direct computation and the fact that 

\x+n-M 



E A + V \x+n-M _ W 
, w w-(l + v)' 
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There are four cases to prove for n > M + 1 and one case for n < M. 

When n < M, direct computation shows the following if we require | — | < 1 then ^2 x>M _ n ^t+u- - - 



w—(l+v) 

n>t n>t 1 I a I a dw (w-a) n ' j e tw w 

< ®n-k >~ (2m)2 d ™ dv w ^ + l _ ^ n _ k+1 Sv+l) w _ (l + ^ 

and if we evaluate the only pole at w = 1 + we have 



then let z — v + 1 — a 



2m J ra _ 1 (v + 1 - ay~ k+1 



2m L ^- fc+1 ^' 



When n > M + 1, 



The cases for (i) j < M, k < M , (ii) j > M + 1, k > M + 1, (Hi) j > M + 1, k > M are 
all similar to the first calculation. 

(i) j < M,k<M : 

Since we have ^^-(a;) = 0ifa;<M — n — 1, 



x>Af-7V 



and after summing up ^2 x> M- n (^f) x+n M = w -(i+z) ^ requiring < 1 and evaluating 
the only pole at w; = z + 1 then 

1 I A I A {* + 1 ~ Q ) M ~ j 1 f ^ 

i^r 7r Q _i /r „ ^ - « (v + 1 - «) M fe+1 



and now we evaluate the pole at z = v 



j-k+ 



2vi Jr^! (v + l-a) 

(ii) j>M + l,k>M + l: 
Similarly to previous proofs we have 

5 The pole at w — does not exist anymore after we summing up x. 
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_n,t 

x& x>M-N 



n-k\ X ) 



since ^>^Ax) = if x < M — n — 1, 

and after summing up ^ x>M _ n ('4r) :E+n_M = w-(i+i>) by requiring < 1 and evaluating 
the only pole at w = v + 1 then 

27TZ J ro ^- fc + X J ' 

(iii) j < M, k > M + 1 : Since we have V^Ax) = 0ifx<M-ra-l, 



xeZ x>M-N 



and after summing up Y^ x >M-n(w') X+n M = w _^ +v ^ by requiring |^| < 1 and evaluating 
the only pole at w = v + 1 then 

' 'v + 1- a) M ~ j = 0. 



2m J ro 



(zv) j>M + l,fc<Mis different. 
Now we require 



then after summing up J2 x>M _ n (^f) x+n ~ M = w _™ l+z ^ an d substituting the pole w = z + 1 
we have 



/* /* ^ — J? 



then it's not hard to see that0 



1 z n ~ j 1 ty n ~^ 

<fc — q(z,v) = dw rrq(w,v) 

z-vz n ~ M y 1 w-vw n ~ M K 1 

if we do the substitution z = —w — A. 
6 Recall 

~ z + v + A(v + A) M -i 

has poles at —A and — A — v, where A is a positive constant such that —A — v is not a pole inside the contour 
of Tq v in the z plane. 
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Therefore 

If 1 z n ~ j 1 / 1 w n ~ j 

® dz rjq(z,v) = (b dw ^q(w,v) 

2m J Tov z-vz n - MyK ' ; 2m J T _ A _ A _ v w-vw n - MyK ' ; 

1 / 1 z n ~ j 

— z — ■ 4> dz rrg(z,v). 

Am J Tn A A z-vz n ~ MijK 1 

u 1 0,v, — A, — A — v 

Now the integrand is of order at most O(j^) as z — > oo since j — M > 1 and therefore the 
integral in z plane is 0. So, 

>=o. 

□ 
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3.3 Construction of Kernel 



Now we are ready to write out our kernel and again for simplicity, we use r pi) .„ )Pi to represent 
a positive oriented contour enclosed the points: p±, ...,pi. 

From previous section, we know 

"2 

#((ni,fi),a;i; (n 2 , t 2 ), x 2 ) = -^MO**))^ x 2 ) + £ (si)*™* (s 2 ). 

k=i 

For ni and n 2 > M + 1, the kernel can be rewritten as 

M 



tf((m, tx), x i; (n 2 , t 2 ), x 2 ) = -^(("i.*!).^^))^!, x 2 ) + ]T 

k=i 

k=M+l 

-(t> {niMUn2M \x 1 ,x 2 ) + ^((m, ix), an; (n 2 , t 2 ), x 2 ) + AT 2 ((ni, tx), a*; (n 2 , t 2 ), x 2 ). 



(25) 



where 



((ni ' tl) ' ( " 2 '' 2)) (xi,x 2 ) = — / 

2m Jt oa 



dw e ^- t2)w l s 



"2 



and 



Kx((ni,tx),xi;(n 2 ,t 2 ),x 2 ) = ^ C-^OC-,^) 

j=M+l 
M 

X 2 ((nx,ix),xx;(n 2 ,U^2) = E^^(^i)^-,(^)- 



Theorem 3.3.1: 



We can also simplify the expression of 0, fTx and i^ 2 which we will denote the simplified 
version of them as <j),K\ and K 2 , respectively. As a result we have the following 

(n 2 ,t 2 ),x 2 ) = -^^^(an, x 2 )+ (26) 
i^i((ni,ii),a;i; (n 2 ,t 2 ),a; 2 ) + K 2 ((n 1 ,t 1 ),x 1 ; (n 2 ,t 2 ),x 2 ) 

where 

Xa) = J_ i ^ eftl ' t2)W h (27) 
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\f f (yj _ \\ni—M e tiw yX2+U2~M ^ 

K 1 ((n 1 ,t 1 ),x 1 ;(n 2 ,t 2 ),x 2 ) = — & dv & dw— — - — — — -- — -,and 



K 2 ((n 1 ,t 1 ),x 1 ; (n 2 ,t w ),x 2 ) 



(27T«) 2 J Tl Jr w x+ni-M+l £ t 2 v ( v _ 1 )n 2 -JVf w _ ^ 

(28) 



1 I dv I dz I dwe hw (w -l) ni ' M Z ^+n 2 -M 2(2-1) + A 1 



(29) 



(27Ti) 3 J Ta J Tlv J Fo W w xi+ni-M ( z _ iy l2 -M e t 2 z ( y _ z _}_ + v - \+ Az~V 

(z-l + A) M v-l + A 

(v-l + A) M (z - 1 + A)(w - a) - (v - 1 + A)(v - a) ' 

Here we require \v\ < \w — 1| for K 1 and \v — a + 1| < \w — a\ for K 2 . 
Proof: 

The proof is by summing up geometric series. 

712 

K l ((n l ,t l ),x 1 ;(n 2 ,t 2 ),x 2 ) = £ ^ j (x 1 )^ j (x 2 ). 

j=M+l 

Since 0((" 1 '* 1 )'( n2 '* 2 ))(x 1 , x 2 ) = if k > n 2 + 1, we can extend the sum from n 2 to oo. So we 
have 

K 1 d dv I dwetlW ( w ~ l ) ni ( X + v) X2+m ~ M w fc 

1 ~~ (9>rri\1 T„ ri in ,„xi+rai-M ^,n 2 +l p t 2 {v+l) 2-~t 



(27ri) 2 / ro J roi w w Xl+ni ~ M v n 2 +i e t 2 (v+i) (u> - l) fc 



1 / / dwe tlW (w -l) ni (l + v) X2+n2 ~ M v M+1 
• dv A 



(2vn) 2 J ro 7r ,i w w xi+ni-M v n 2 +i e t 2 (v+i) ( w - \)M( W _ i _ v y 

Here we require |t>| < |u> — 1| which can be achieved if we make v small enough. 
After change of variable v — > v — 1 we have our 

^ 1 I dv I dwe tlW (w - l) ni ^2+™2- M 



(30) 



(27U) 2 / ri / ro A W wxi+ni-M ( V - \)n 2 +l e t 2 v ( w _ ( w _ y j 

Now we proceed to K 2 . 
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M 

k 2 ((n u tl ), Xl ; (n 2 , t w ),x 2 ) = £ Kf-MWn-f 0*0 

cW lU, (u> - l)"i- M 



2vn / r 



1 O.CK 



X 



dv j J (l + z )^+™2- M ( w -a) M 2,2 + A 1 



(27ri) 2 J Tal v - a + 1 jf^ z n 2 -M e t 2 ( z +i) ( v + i _ a )Af 2 + v + 

(z + A) M y> (v + A) k (v + l-a) k 
l v + A ) M j~i ( w ~ a ) k ( z + A f 

After summing up the series and make z , w big enough and v close to a — 1 enough such 
that 

\v + l-a\\v + A\<\z + A\\w-a\. (31) 

we have the following 



dwe tiw (w - l)«i- M 

(27ri) 3 J^ Q a W w xi+ni-M 
dv I (l + z )X2+n 2 -M 2Z + A 1 



K 2 ((ni,t 1 ),x 1 ;(n 2 ,t w ),x 2 ) = jf — - .^Jj^ x 



f - a + 1 / ro u 2 n 2 -M e t 2 ( 2 +l) ( v + 1 _ qi)M ■j + u + ^^.y 

(z + A) M (u + A)(u + l-a) 

(u + A) M {z + A)(w - a) - (v + A)(v + 1 - a) ' 



(32) 



If we compute the simple pole at w = a + - (frf^j ^ we have (z + A) (w — a) = (v + A) (v + 
1 — a), then the residue is zero since the contour integral on the v plane is zero. 

Now for convenience which we will see later, we change the order of integration along with 
the change of variables v — > v — 1 and z — > z — 1 then we have 

K 2 ({n 1 ,t 1 ),xi,{n 2 ,t w ),x 2 ) = ,J~ / x 

(2-niy J Fa v -a 

dwe tlW (w -l) ni ~ M z ^+n 2 -M ( W - a )M 2(z-l)+A 1 



dz 



Flv J To W w xi+m-M ( z _ iy 2 -M e t2Z ( v _ a )M z _ i +v _ i + a z _ v 

(z-l + A) M {v-l + A){v-a) 

(v-l + A) M (z - 1 + A)(w - a) - (v - a)(v - 1 + A) ' 

(33) 



□ 
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4 Asymptotics 



In the next three sections, we will sketch the proof of convergence of the rescaled process Xt 
in each region in the sense of finite dimensional distributions. We will follow the procedures 
described in [3]. 



4.1 Slow Particles Region 

Theorem 4.1.1 

Here we will show 

lim X t (u) = DBM(— ln(a(v))) 7 v e (0, (1 - a) 2 ) (34) 

t— >-oo 

in the sense of finite dimensional distributions. It should be noted here our results is similar 
to 0. 

To prove it, we prove a more general version of Theorem 3.4.1. 
Proposition 4.1.2. 

We let 7r(0) be a function from K -)■ R such that < 1. Define U = (n(6i) + t )T and 

u' 2 -2a+2 n 
a(2-a) P ' 



m = (tt(^) - 6 t )T + M then for < n < (1 - a 2 )t which means < tt(0) < a V 2 2 _ a t 2 g, we 



have 



x n(gT )(t(e,T)) - (act-2=^) 

X T {9) = { ' n V ;7 ^ ->• DBM(r(0)) as T -> oo (35) 

—ayT 



with a 2 := a(vr + 0) - and r(0) = - ln(a). 

Proof: 

We need to show the kernel 

^((Tii.tO^i; (n2,t 2 ),a; 2 ) = -0 (ni '* l} ' (n2 '' 2) (xi, x 2 ) + + K 2 ]{(n 1 ,t 1 ),x 1 ; (n 2 ,t 2 ),x 2 ) 
converges to K DBM after rescaling. Here we have 



xt = orti - rij/(l - a) - £i<JiVT, (36) 
u . = vr (^.) + 0h (37) 

0i = 7r(0i) - e h (38) 
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of = aui — aaj/ (1 — a) 2 . 



(39) 



*i = ihT, (40) 

n^M + a/T. (41) 

Since we require t\ > t 2 if n 2 > rii, we have a 2 — ai > and — w 2 > with one of them 
being strict inequality. 

There are three parts to show. 

(6 - 60-2M) 2 



cxp 



lim Cd> 

T-s.00 



2(1 - ot/ol)) 



lim Ci^! = 0, 

r r W M 1 p w 2 /2+^i 

Jim CK 2 = f dV & dW- 

J\V\=R J-i 



t^oo J ]vl=R J_ L+lR V M W - V e y2 Wi) a +^M 

with L > R and 

r; = — 

Q,a;2— a;i+n2— ni 



„ W1-U2 a 2 - ai 
a (1 — a)^a 



<£(«) 
02(a) 



ct 

1 



(42) 



c = ^— ^ ^Vr. (43) 



Analysis of 0: 

First, we define the following functions 

a 

go(w) = (ui - u 2 )(w - a\n(w)) + (ai - a 2 )(ln(u; - 1) + hi(u>)), 

1 — a 

9i(w) = ^ 1 a 1 -^a 2 )Hw), (44) 
g 2 (w) = - ln(w). 

In order to do steepest descent, we must calculate their derivatives and we found 



#d(a) = °> 

Ui 
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Using w — a = ae iy — a = a(iy + \{iy) 2 + ■■■) their Taylor expansion at w = a are given as 

v 2 

g (w) = g (a) - a\) + 0(y 3 ), 

9i(w) = gi(a) + iy(£i<Ti - &a 2 ) + 0(y 2 ), ^ 
g 2 (w) = - \n(a) + 0(y). 

Let T = {ae iy , y G (— it, it]} and use e x ~ 1 + x + 0(x 2 ) for x « 1 then for some 5 > we 
have 



<f>= <f dwe T9o{w)+VT9l{w)+92{w) 

= I ( j iWe T9o(u>)+VTg 1 (w)+g 2 (w) + / ^ we Tgo(w)+VT gi (w)+g 2 (w) 

J\w-a\<5 Jr \{\w-a\<8} 

= I dwe Tgo(w)+V?Tgi{w)+92{w) + e T 9o( a )+^9i{^)+92(a) I dwe~>* T , 

J\w-a\<5 Jr \{\w-a\<6} 

After multipling by C the second term has an error of order 0(e~ MT ) and fi ~ 5 2 . After we 
drop the error term, we have 

= I rfwe T( 90 («)-|y 2 ( CT r^))+v^(9i(«)+^(6^i-6^))-ln(a)+TO(^)+v / TOfe 2 )+O( ?/ ) 

ir n{|w)-a|<5} 

= e Tg {a)+VT gi (a)-ln( a ) I ^-T^yHaf-a^+iyVT^ai-^) x e TO(y 3 )+VTO(y 2 )+0(y) 

Jr n{\w-a\<8} 

= e Tg (a)+VT 9l ( a )-\n(a) I ^^-T^^l-aD+iyVTi^-^) 

ir n{|io-a|<<5} 

+ e Tg ( a )+VT gi (a)-ln(a) I ^Q-T^^-^+iyy/T^un-^) ( e TO(y 3 )+VTO(y 2 )+0(y) _ 

Jr n{\w-a\<5} 

(48) 

Now we use \e x — 1| < el x l|x|, then we have the second term of the above equation becomes 

cnst . | I ^ e -T|y 2 (af-a2)+^v^(aa 1 -6^2)( e TO(y 3 )+v / TO(j/ 2 )+0(s,) _ -Q | 

J r n{\w -a\<6} 
r 5 

< cnst / ^| e - T ^ 2 ( a i- a 2)+^«i CT i-?2^)( e TO(s/ 3 )+VTO(j/ 2 )+o(s/) _ ^| ( 49 ) 

Jo 

< cnst / ^|e- T ^ 2 ^-^)+^^^-^- 2 )|| e o(^+^ 2 +J /) 0(2/ 3 + ^/f y 2 + 

Jo 

We do the change of variable p = yv^, then the above becomes 
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cnst 



rVTS 

/ e - cp2+ -0(p 3 + ...)dp. 
Jo 



After mutipling by y/T, the estimate of error is 0(1/ \/T) which goes to as T goes to oo. 
So now the last step is to calculate the first term in equation (4.1.22). 

e Tg (a)+VT gi (a)-ln(a) I ^^-T^iaf-cr^+iyVT^cn-^) _ 

Jr n{\w-a\<5} 

If we let l/D = %T + l\ = e T ^) + Vr 9l{ a) , w = as*, , z/Vf = y , and extend 

5 — > oo for w G Tg = {w G r , |w — a\ < 5, w — a = re ±Mr / 2 } for some r then 

D /r n{|w-a|<<J} 

= l/(DVf) [ dwe-^ z2 ^- a ^ +iz ^ ai -^ + l/(DVf)0(e~ 5 ' 2T ) 
Jr 

= 1/C [ dwe-**rt-^ ) **K lOl -t an) + l/(DVf)0(e- g ' T ). 
Jr 

Thus when we mutiply the above equation by C we have 

1 (Sl-C2^2/^l) 2 
1 — T 2,„2* 



as T — >• oo. 
Analysis of A^ : 
Similar to previous proof, 
we rewrite 



Ki({n 1 ,t 1 ),x 1 ;{n 2 ,t 2 ),x 2 ) = f dw <p dv 

Jr, Jr n 



f/o,iW+v^/i,iH i 
r Fl j Fq e Tfo,2(v)+Vrfi,2(v) w(w - v) 
where / 0)i = «i(u> - a ln(iu)) + aj(ln(w - 1) + ^ ln(iu)), and / M = &<7j ln(iu). 

The derivatives of / ,j an d are 

,, , x _ (w - a)(ui(l - a)(w - 1) + di) 
u>(-u7 — — a) 

ell I \ OLUi .1 OL . 

W = % > o. 

cr 

/l» = — • 
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(50) 



C-(j)^ —== e w-H'V (51) 



(52) 



There are two critical points for Jq^w), namely 

W- = a and w + j = 1 — - — -. (53) 

Ui(l-a) 

By using the inequality rii < (1 — a) 2 t, ti = u{F and rii = M + diT, we have 



M + aiti/u)i < (1 - a)\ 
=>- di/ui < (1 - a) 2 < 1 

1 - > a > 

1 — a 

=4> < w_ < < 1. (54) 

We choose r = {ae^, e [— n, n)} and I\ = {t> = 1 — re 1 ^, r = 1 — w +2 , ^ G [— 7r, 7r)}. 
Observing \w — 1| 2 > (1 — a) 2 > ai/ui we have 

d$t,(fo t i(w)) au\ sin< 



|2 



(|w-l| -ai/ui), 



o?0 |w — 1| 

which is decreasing away from w = a and from \v | 2 > |u> +2 | > |-u7 +2 | 2 a: we have 

d$t(-fo i2 (v)) ru 2 simp 2 _ a 2 „ 

# " " M 2 m a[ u 2 {l-a )h 

which is decreasing away from v = w +2 . 
So now we have 

^i((ni,ti),Xi; (n 2 ,t 2 ),x 2 ) 
j/ciW+Vf/i,! 



= (j) dw (j) 



dv 



Tl J ro e Tf ,2{v)Wrh,2(v) w(w - v) 

If dw f exp(T%{/o, 1 ( W -)-/o, 1 (a)} + ...) 1 
C/r/ 7r exp(TR{/ 0>2 («; +2 ) - /0.2(a)} + ...)«;(«;-«) 

' ' du> <f dve~ ST+ - ' 



C 7ri 7r «;(«;-«) 
The leading term has an order of e~ 5T for some 5 > thus 

Ctfi ->■ as T ->■ 00. 

Analysis of K 2 : 

Again after substitution we have 
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K 2 ((ni,ti),xi, (n 2 ,t w ),x 2 ) = j> dv (h dz (b 



dw (w - a) M 
w (v — a) M 

e Tjb,i(«)+Vr/i 1 iM 2 (z - 1) + A 1 
e T/ , 2 (*)+>/57i,2« z-l + y-l + Az-v' 

{z-l + A) M v-l + A 

{v-l + A) M (z — l + A)(w — a) — (v — a)(v — l + A)' 



(55) 



Now we study the above formula in two cases: 
Pole at z — 1. 

The steepest descent path are 
{u = a + fe^,^G [0,2tt)}, 

{w = (ct — e)e i( ^, G [0, 27r)} since w = a is a root for the integrand, and also a saddle point, 

= l _ re in , r = 1 - w%, Q G [0, 2tt)}. 
For e small enough the leading contribution is ^e~ 5T for some 5 > 0. So as T — > oo, the 
contribution from the pole at z — 1 goes to zero. 

Pole at 2 = t). 

Integrating the simple pole at z = v we will have 



K 2 ((n 1 ,t 1 ),x 1 ;(n 2 ,t w ),x 2 ) ^ S dv (b 

JT a J Tq 

T/ ,i(i«)+v / T/i,i (to) i 



diu (u> — a) M w M 
~w~ (v-a) M v M 



e Tfo,2(v)+VTf 1>2 (v) W - V ' 



(56) 

Now we take the following steepest descent path 

T a = {v = a + R/Vfe 1 ^^ G [0,2tt)} 

and 

r = {w = (a - L/Vfy 41 , G [0, 2tt), L > R}. 

Using Taylor series of 

/o» = /o» + + 0((w - a) 3 ), 

/i» = / M (a) + ^ W ~ a) + 0(( W - a) 2 ), ^ 
(w(w - a)) M = a M (w - a) M + 0((w - a) M+1 ) 
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and the change of variables v — a — aVa 1 l T 1 I 2 ) and w — a = aWa 1 1 T 1 / 2 we have 

Tf 0A (W) = T/ ,i(a) + ^ + 0(1/VT), 
Vff ltl (W) = Vff 1A (a) + WCi + 0(1/Vf), 

Tf , 2 (V) = Tf , 2 (a) + ^ + 0(1/Vf), (5g) 

Vff lt2 (V) = / 1)2 (a) + + o(i/Vf), 

(w(w - a)) M _ + 

~ V M l + 0(^)' 

Applying the above substitutions, we have 

. . f f dw(w-a) M w M e T ^ (w ^^^ 1 
K 2 {(n 1 ,t 1 ),x 1 ;{n2,t w ),x 2 ) = f dv (b — 

JT a Jr W \ V ~ OL) M V M e r /o,2(u)+VT/i i2 (u) w — V 



r W 2 /2+W^ W M i 

dV i dW ., 2 , . - — ^—r— — - + 0(1/Vf). (59) 

Here L > R. 



So now we have 

r r e w 2 /2+w^ w m 1 

C ■ K 2 ((n 1: tl ), Xl ; (na, t w ),x 2 ) -> ^ £^ ^ ^/^m 

as T — >■ oo. 
□ 
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4.2 Linearly Decreasing Region 



Theorem 4.2.1 

a 2 -2a+2 



For n e ((1 - a) 2 )t,t) or " a(2 2a + 2 < n(9), The rescaled process 



converges to the Airy 2 



Xr(0) = -H-C-JV^) (60) 



hm X T (r) = S v A 2 (r/S h ). (61) 

T-s-oo 



Proof: 

To prove this theorem, we have to show after rescaling our kernel K = 4> + K\ + K 2 — > K.Airy 2 - 
Similar to previous theorem, there are three parts to show 

[( i2 _4)-(4 L -4>] 2 

°v St St 

h h 



lim Cd> = , 






/ 4 "<f ' 




lim CK X = 


/ dV J 


71 


lim Cir 2 = 

T^oo 


o, 





T 2 Tl 



i {r2>ri}) 

v 3 /3+r 2 5 h y 2 -v(5 2 / M - K 2 T |/ K0 )/ K j/ 3 1 (62) 



e [/3/3+ri5 /t [/ 2 -C/(si/ M -K?r 1 2 / K0 )/4 /3 [/ - V 



and the conjugation constant C = e to-^ J + i// 2 ' 1 ' J /tnWT 1 / 3 . 
First of all, our rescalings are 

t- = T[tt(# - r,T- 1/3 ) + 9 - nT- 1/3 } ~uT- r t (l + tt')T 2/3 + Vr 2 T 2/3 , (63) 



ni -M = n(9- r.T- 1 / 3 ) + r, + T 2 / 3 (l - tt') ~ oT + r,(l - tt')T 2 / 3 + Vr 2 T 2 / 3 , (64) 
Using Vl + Ax = 1 + ^ - |A 2 x + • • • , 



Xi = ti- 2y/ti( ni - M) - s,(T 1/3 ) 



= u(l - 2v^)T - r,T 2 / 3 [(l + tt>* + j— ^ 
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T ^ mw - - jj-) + " 1 ^;^->' i - «cm ( 6 5, 

Analysis of 0: 

If we write = ^ e 9i(« , ) T2/3 +32(>)T' 1 / 3 j^ then assume r 2 > Ti for the rest of calculations 

1 - 7T' 

ffi(w) = (ri - r 2 )[-(l + tt')w + (1 - tt') ln((w - l)/iu) + [(1 + tt')w* + ]lnw] 

1 — w* 

= »(»•) - i.n - r>)\(^ - w ,*l~_l, f )(u - »-) 2 + 0({w - to*)) 3 , («> 

= 9 iK) - (t, - r 2 )K,(m - ti>*) 2 + 0((w - w")f. 



+ (si - s 2 ) lnw 

= #2 ™ - b T l _T 2 —\{W-W 

2(1 — 2w* u w* 

+ 0((w-w*f) 

= 92(w*) - (^(7? - rl) - S -^)(w - w*) + 0((w - w*) 2 )- 

(67) 



So we have 



0=1 ^ e < ? i(»)T 2 / 3 +< ?2 (»)Ti/3 
Jt w 



2 

3 



e T 2 / 3 9l («,*)+TV3 92(lu »)_ T 2/3 (ri _ T2)Kl(w _ w . ) 2_ T l/3Q (r 2_ T 2 ) _£J^ )(w _ w , H ^^^ 

r 

(68) 



Let (w — w*)T 1 / z k 1 J^ = iy and using J R e ax2 / 2 + lJx dx = \J^fe j2 /( 2a \ then it becomes 

2tt y R 



^2 T 2 
rf s 2 _ T 2 w s l _ T l ni2 
^ 
h 



= e TV3 gi( ^ )+r i/a g2(w . )T _ 1/3 -1/3 ,-1 / 1 ^ e iT^T) 0(T -l/3 } 

V 47r( t-it) 

where = k^ 1 ^ 3 and = w*k^ 3 . 



(69) 
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Analysis of K\\ 



f (w) = wu + a\n((w - l)/w) - [u(l - 2y/aju) + a]\nw (70) 

/;W=M( l + ^---l^) (71) 
w — 1 w w 

One can check at w* = 1 — s/u is a double pole such that /q(w*) = fo(w*) = 0. 

/oH = /oK) - ^ " - ^*) 3 + 0((«; - ^*) 4 ). 

OW*{L — W*) 



fi,i( w ) = + + (1 - *r>i Hi™ - l )l w ) + + + ^3^:] (72) 
The Taylor series of the above functions are of the following: 

/i» = hM) - - w ^Zly )(™ - -*) 2 + 0((w - w*) 3 ) (73) 



r I \ rl ////>W i / W— 1\\ rl litres! 1 \ ((1 ~~ U ) 7r ' ~ (1 + Z/ )) 2 m 2 i 

/ 2 » = [-tt"(0)(™ + ln( — )) - [-n"(9)(a - — ) + ^ _ K + ^ ^ 

(74) 



/,(«,) = feW-K / ^ Lffi ^ 1 ^ ^--^]( W -^)+Q(( W -^) 2 ) (75) 
If we let k = *~7T~~ — 1\ and «i = iU^f^ rrr^r^), 

u u>*(l— u>*) 1 2 v ui* mj*)^ / ' 

then 

1 

/to(w — w 

/i,iW = /i,iK) - - u>*) 2 + - W *) 3 ), 



/oH = /oK) - - ^*) 3 + 0((w - ^*) 4 ), 



(76) 

/ 2 » = / 2 >*) - (^r 2 - - «,•) + 0((w - w*) 2 )- 

Kq W 
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With the rescalings, and substitute the above functions, we have 



p - T 3^^n^( W -^) 3 - T2/V 2 K l(^- W *) 2 - Tl/a (^- S 2/«'*)(«'-^) 1 

K x = B- 1 & dv & dw 



B 



T 3^n"-^*) ( ^ w * )3 ~ r2/3TlKl ^~ w * )2 ~ T1/3(!! lr~ Sl/a)(t '~"' ,) w ( w ~ ^ ' 

(77) 

1 _ e T 2 / 3 /i i2 (»*)+T 1/3 /2,2(™*) 

Let (w — w*)(/t T) 1 / 3 = w 2 and (v — w*)(/t T) 1 / 3 = w 1: then the above becomes 

r r f> w^+w^T 2 /s h -(s 2 /s v -T^/sl)w2 i 
l/C / dw 2 / d Wl — 2/c - - — „ + 0(T- 1/3 ). (78) 

Analysis of K 2 : 

After inserting the rescaling terms, we have 



r r r e -r^ I ^(i»-^) 3 -T 2 /v 2K1 ( w - w *) 2 -TV3 ( ^_ S2/w . )(u ,_ w * ) 1 

A 2 = -B" 1 6 dv (h dz <£ dw — z r 

JV 1 Jl,v Jr ^T ^^^ z-w^-TV^K^z-w^-T^C^-si/^iz-w*) W{Z - V) 



{w-a) M 2(s-l) +4 {z-l + A) M v-l + A 



(v - a) M z - 1 + v - 1 + A (v-l + A) M (z - 1 + A)(w - a) - (v - 1 + A)(v - a) 

(79) 



B 



1 _ e T^ a h,2i w )+ T ' ' J /2,2(« <*) 



e T 2 /3/ lil ( u ,*)+T 1 /3/ 21 ( tu *) • 

The leading comes from the T -1 / 3 — neighboorhood of w*, so we let z = w* +Wi(T 'kq)^ 1 ^ , w = 
w* + w 2 (Tk )~ 1/3 , then 



K 2 = T^K-^B- 1 I dv I dw 2 I dw 1 



e w'i/3+w 2 T 2 /S h -(s 2 /S v -Tl/Sl)w2 

,r\' I! 4-ir-- S', — ( ci S - S'r I / / ' 



Qwf/S+wln/Sh-isi/Sy-rl/SDwi 



w*-a) M 2(w*-l)+A (w*-l + A) M v-l + A 



(v - a) M w*-l + v-l + A (v-l + A) M (w*-l + A)(w* - a) - (v - 1 + A)(v - a) ' 

(80) 

After multiplying by the conjugating constant, it goes to zero as T goes to infinity. 
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4.3 Transition Process 



Theorem 4.3.1 

2-2Q+Q 2 



When n ~ (1 — a) 2 t, n(9) = 2 a ^ 2 ^) Q with 9 fixed and we let the rescaled process be 

x n {t) - (t - 2y/t(n-M)) 



X T(r)- _ Tl/3 
We would like to show 

lim X t {t) = S v A DBM ^ 2 (r/S h ), (82) 

where Sh = kJ 1 ^ 3 and S v = w*k^ 3 . 
The process A DBM _> 2 also appeared in [2]. 
Definition: 

The Aobm^2 process is an m— point distributions at T\ < r 2 < ... < r m such that 

H^ =1 A DBM ^ 2 (T k ) < s k ) = det(I - X»KA D BAt->2Xa)v>(n,...,T m )xK) (83) 
where Xs(r k ,x) = I x>Sk and 

= - vm^) 1{r2>ri} + A A r/ "'v^-- - --v/',,, „- 2 - 

(84) 

Here we have 71 : e~ 27r / 3 oo — >■ e 2? ™/ 3 oo and 72 : e^^oo — > e~ m ^oo provided that they both 
pass on the the left of and they don't meet each other. 

Proof: 

To prove the convergence, we have to show after rescaling our kernel, K = <p + K\ + K 2 — > 
Similar to previous theorem, there are three parts to show 



[( £2_4)-(iL_4)] 2 

a/Qi„„2 TO /c _/ „/Q _^2/q2 



lim CK\= / c?w 2 / rf^i 



e TO|/3+w|r 2 /Sh-(s2/S^-r|/S^)it)2 ]_ W M 



lim Ci^2 = <p du dw 2 / (iwi — - — 5 — 7- — : — — — — -. —. - — 

T^oo J To J r2 J 7i ^wf/S+wfTx/SH-isx/Sv-Tl/SDwx ^ - U )(w 2 - u) U 

HprP r = e Q(t2 - fl) (Q-l)" 2 ' ni T l/3 

lid C v_ rv x 2 — x l + n 2 — n l 
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From previous section, we have 

ti = T[n(9 - nT- 1 / 3 ) + 9- nT- 1 / 3 ] ~ uT - n {l + tt')T 2 / 3 + Vr 2 T 2 / 3 , 
m - M = n(9 - TiT- 1 ' 3 ) + n + T 2 / 3 (l - tt') ~ aT + nil - tx')T 2 ' 3 + Vr 2 T 2 / 3 , 

1 - 7r' , 



Xi = u{l - 2^Ju)T - r,T 2/3 [(l + n')w* + 

, r 2 T l/3 ((l-^K-(l + ^)) 2 _ fT l/3x 

+ * 4u(l-«;*)3 J 

= D T - DmT 2 ' 3 + (L> 2 r 2 - Si )T 1 ' 3 . 

Also 



1 — w* 



w* = 1 — = 1 — (1 — a) = a. 
If we write = ^ ^e^ 1 ^> T ' 2/3 +^^) rl/3 [Jhen 



1 - 7T 7 

= (ri - r 2 )[-(l + 7r')w + (1 - tt') ln((w - l)/w) + [(1 + 7r')a + ] law] 

1 — a 

1 1 + 7r' 7r' — 1 

= - (n - r 2 )-( - a) 2 + 0((w - a)) 3 , 

2 a a{± — a) z 

= gi(a) - (n - t 2 )kx(w - a) 2 + - a)) 3 . 



/ \ A H f n\ / , r l »,„ w 1 v. ((1 - VW - (1 + Z')) 2 , 2 

2 to 2 1 — a 4tt(l — a) * 

+ (si - s 2 ) \nw 

= »(<*) - (^(t? - r 2 2 ) - - a) + 0((w - a) 2 ). 

Oi 



So we have 



^ e9l HT 2 /3 +g2(M , )T l/3 
l W 



e 
r 



T 2 /3 9l ( Q ,)+rV3 92 ( Q ,)_T 2 /3(r 1 -r 2 )Ki(u>-a) 2 -T 1 / 3 (^(r2-r|)-^^)(^-a)+T^ 



7 When tt(9) = 2 ~ 2 "+" (9 

v / a(2 — a) 
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Let (w — a)T 1 / 3 Kg^ 3 = iy and using j R e ax2 l 2+tJx dx = y^fe j2 /( 2a )^ then it becomes 

27T Jr 



■ "2 _ T 2 \_ t s l _ T l n,2 



= e ^ l(a)+ TV 3g2 ( a ) r i/3 K -l/3 a -i / 1 e - 4^-^ +0(T -l/3 } 

V 47r 'sh ?J 

where = K^rif/ 3 and S*„ = cxk^/ 3 . 

Now we are going to show the rescaled kernel of K\ converges to 

Iff f e wl+wlT 2 /S h -{s 2 /S v -Tl/Sl)w2 l 

du / dw2 / dw\- 



e T 2 /3 /l l(a)+T i/3 /2 i(q) • 

Let (iu — ct )(Tk ) 1/3 = w 2 and (f — a)(T/t ) 1 / 3 = then the above becomes 

If f e ™l+™ll-2/S h -(s 2 /S v -T%/Sl)w2 i 



C- 1 = B-^T-Wa- 1 



Similarly, the rescaled kernel of K 2 converges to 

Iff f e w^+w? 2 T 2 /S h -(s 2 /S v -T^/Sl)w 2 i Wi 

du / dw2 / dwi — ~ y— — - — — — 5—^7 — — -( — ) 



(27ri) 3 Jr o 7 72 7 7l e ^?+-ri«)f/5 h -( S i/s l ,-r 1 2 /52) Wl ( Wl _ M )( W2 — u) u 
as T — >■ 00. 



(91) 



(27ri) 2 / ro 7 72 J 71 1 e t^+n«;?/s h -( a i/s„-T?/s2)«'i Wl _ W2 ■ 
We first write 

f f e TfoM+T 2 / 3 f 1 .,( W )+T 1 /- i f 2 ^(w) + ... ]_ 

Kl = 1 ^ 1 dW e ^W + ^A. < W + T^, <W+ - w{w - v) ^ 
Similarly to previous section, the leading term is 

f r e -T3^j^( W -a)3-T 2 /»r 2Kl ( W -a) 2 -TV3 ( ^M_ S2/a)(w _ a) 

i^i = B' 1 i dv <p dw — 2 r, (93) 

Jy 1 Jr e -T^^{v-a^-TV^T 1 K 1 (v-ay-T^('^- Sl /a)(v-a) W{W - V) 

D _i p r 2 /3 /l2(a)+T i/3 /2j2(a) 
X5 = 



/ ^2 / ^l 1 ! 2lo , la 27^ ~ +0(T- 1/3 ), (94) 
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5 Appendix 



5.1 Dyson's Brownian Motion 

In this section we define what a Dyson's Brownian Motion (DBM) is [10] .We introduced it 
from a random matrix point of view [10]. Suppose, we have a matrix H^xn with eigenvalues 
Xj, 1 < j < N, then the degree of freedom of the matrix is d = N + N(N — l)/3/2 since the 
matrix H is determined by its independent real parameters which are H ii} 1 < N; H^\ 1 < 
i < j < N, < A < /3 — 1. Now we rename other parameters with with [i runs from 1 
to d and use it to replace X,i, and j0 Suppose at time t, they have values Hi, ...,Hd and 
Hi + 5 Hi, Hd + SHj at time t + 5t. Brownian motion is defined by saying each infinitesimal 
5H^ is a random variable such that 

<5H„> f = -H^St (95) 

< (5H,) 2 >f = g,kT5t (96) 
for some given constants /, k and T. Here 

^ = ^ = 1 + ^, (97) 

All other averages have a higher order in 5t. The Smoluchowski equation which corresponds 
to the above requirements is 




The solution P(H\, ...,Hd,t) to the above equation is the time dependent joint probability 
density of H^. If H = Hi at time t = T\, H = H 2 at time t = T\,f = 1/2, 1/(3 = 2 = kt and 
H are M x M hermitian matrices (N = M) then we have the following formula[10j. 

Definition 5.1.1. 

P(n, Hi; r 2 , H 2 ) = — _ ^_ n)))M ^(- ^Se-L-n)) } ( " } 
The joint distribution of the largest eigenvalue of the stationary DBM process are given by 

P(nr =1 DBM(r fc ) < Sfc ) = det(I - XsK^Xs)^,.,^}^) (100) 
s H fl are uncoupled and each is subjected to a fixed simple harmonic force. 
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Here the kernel is given by 

/ (x 2 -x 1 e-^- T ^) 2 \ 
exp I — — _ t m-1 

K»™(n,XM,x 2 ) = ylJ^^ + £ ^^PkMPkMe-*" 

VMl - e-2(^-n)) ^ 

(101) 

where = H k (x/V2)n x / 4 2 k ^ 2 (k\) 1 / 2 , and H k (x) is the standard Hermite polynomial 

of degree . -Kdbm is also called extended hermite kernel. 
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5.2 Extended Hermite Kernel, Extended Airy Kernel and DBM 
to 2 Kernel 



In this section, we give an explanation about the double integral representation of the ex- 
tended Hermite kernel of the Dyson's Brownian motion K^bm- 

The extended Hermite kernel is defined as 

K DBM (t,x;s,y) = — — exp{-( — - )}xt,s + ^2tyq k p k (x)p k (y)e- y (102) 



k=0 



where p k (x) are normalized Hermite polynomials, Xt,s — \<s an d q = e l s . p k (x) can be 
written as p k (x) = . 1 H n {x) [?] where 

v / 7r2 n n! 

HJy) = ^=e y2 [ dwe w2 - 2wy w n . (104) 

Here the contour of Y = {w e -L + iR} and 7 = {z = re i27r< ^, G [0, l),r < L} with L,r 
being some fixed positive constants. 

Now we know H n (x) — if n < — 1, we can extend the following sum to —00, that is 

rt— 1 00 
k=0 k=— 00 

l^J 2 ,if ] T dW % dZe 
After summing up the geometric series 

^ z k w — z/q z n ~ l 

k=— 00 



by requiring \w/(zq)\ < 1, we have 



J]q k Pk{x)p k {y)e- y2 = -^— 2 [ dw I dze w2 - 2yw -* 2+2xz 1 ^ 
^ (2tu) 2 J v J 1 w- z/q z r - 

Now we use the change of variable z — > zq then we have 

Tq k Pk(x)p k (y)e- y2 = [ dw I dze ^-^-^ qxz (1Q5) 

^ (27r«) 2 J r J 7 w - z z n 
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If we let w — > —w/y/2 , z — > —z/^/2 x — > —x/\^2, and y — > —y/^/2 we will have 



y2q k p k (x)Pk(y)e y2 = , 1 . <f> dz [ dwe 

( 2m ) J\z\=r J L+iM. 



1 in 11 
w 2 /2+yw-q 2 z 2 /2-zxq 1 w 

W — Z Z n 



Thus, an other way to express Kdbm is 



K DBM (t,x;s,y) - ^ ^l ^ lfo- 



dz / dwe w 2 /2+y W - q 2 Z y2-zx q _ 
\z\=r JL+iR 



(2ni) 2 J 

which is what we had in previous section. 
If we use 

^T^-pM^^)}^ = f>(*)Mt)/e- 2 , 
^ y > y fc=0 

we also an alternative formula 

j _ E - nPfe (,) PfeW ^e- 2 , t<s. 

Similarly we can also rewrite the extended Airy kernel 

v ( \ f Jo°° e^~^Ai( Sl + A)Ai(s 2 + X)d\, n > r 2 , 

^(r 1>ai ;r 2 , a2 )-| _ e -,(n-. 2 ) Ai(si + A)Ai(s2 + A) ^ n < Ta . 

as the following 

rf-r? , , 1 [( S2 -r|)-( 81 -r2)] 2 

e __ (T2S2 _ TlSl) ^ (ri)Si . r2)S2) = _ ___ e J{T2>Ti} 

If f e W%+T 2 W%-(s2-T 2 )w 2 ]_ 



(2vri) 2 J 72 J 71 e u>i+T 1 w 2 -(s 1 -r 2 ) Wl Wl _ W2 

by using 

Wi 



1 /"°° 

= - / e x ^~ w ^d\ 

- w 2 Jo 



^^7t(t 2 ~ Ti) 



e 4 (-2-i) 2 + 12 = / e A(r2 - Tl) Ai(s 2 + A)Ai(si + A)rfA 

</ — oo 
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from, [6J and 



^ J dwe w3/3+aw2+bw = Ai(a 2 - b)e 2a3/3 - ab . (114) 



Now we show how the kernel Kdbm->2 in the proof of previous section can be written in the 
form o£] 



e 4(73-71) if f e w^/3+T2W^-s 2 w 2 W M J 

I + 7^Z7^ / dw 2 / rf^i -375T7-2-iT- -77- — • (115) 



v/47r(r 2 - n ) (2™) 2 7 72 V 71 1 e ^»Hi»i<«; 1 -tt; 2 
The key is to use 

t w ^ = (l/(w 2 -u)-l/(w 1 -u)) (116) 

{w x -u){w 2 -u) w 1 -w 2 K n J n JJ K J 

and S=S = ^Er=o( M / w )" ; b y letting |«| < |uw|. 

Once we use the identity and integrate out the integral on the u plane, cancelation happens 
and we have 

[(&-i>-(&-^)] 2 



K 2 n h — i /• r e wl+wlT2/s v -(s 2 /s v -T^/sl)w2 i 



4vr(f-f) s (2fi) 2 L 2 L 1 e ^»?n/s t -( S i/s,-r 1 2 /sD'«i Wl - W2 



/" /" e ^2+ w 2 r 2/5 , h -(s 2 /5„-r 2 2 /52)u;2 I W M 

i I*' > 3 Jr ^ iy 2 ^ L ^ e wi+w2 ^ s *-^ s ^/sl)^ ( Wl -«)( W2 _ U ) ^ 



1 dfH 



[( _j2 _4)-(^-4)] 2 

ii h 



4vr( 



6 h 6 h 



T£_ _ Tl\° Ji T 2 >Tl 

s h s h ' 



e w^+w^T 2 /S h ~{s 2 /S v -r^/Sl)w2 W M 



(2vn) 2 y z 7 1 e^+^n/^-lBi/^-^/s^K w m Wi _ W2 ■ 



(117) 
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